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MULTIVARIATE SUPOU PROCESSES 

By Ole Eiler Barndorff-Nielsen and Robert Stelzer 

Arhus University and Technische Universitdt Miinchen 

Univariate superpositions of Ornstein-Uhlenbeck-type processes 
(OU), called supOU processes, provide a class of continuous time 
processes capable of exhibiting long memory behavior. This paper in- 
troduces multivariate supOU processes and gives conditions for their 
existence and finiteness of moments. Moreover, the second-order mo- 
ment structure is explicitly calculated, and examples exhibit the pos- 
sibility of long-range dependence. 

Our supOU processes are defined via homogeneous and factoriz- 
able Levy bases. We show that the behavior of supOU processes is 
particularly nice when the mean reversion parameter is restricted to 
normal matrices and especially to strictly negative definite ones. 

For finite variation Levy bases we are able to give conditions for 
supOU processes to have locally bounded cadlag paths of finite vari- 
ation and to show an analogue of the stochastic differential equation 
of OU-type processes, which has been suggested in [2] in the univari- 
ate case. Finally, as an important special case, we introduce positive 
semi-definite supOU processes, and we discuss the relevance of mul- 
tivariate supOU processes in applications. 

1. Introduction. Levy-driven Ornstein-Uhlenbeck-type processes (OU) 
are extensively used in applications as elements in continuous time models 
for observed time series. One area where they are often applied is mathemat- 
ical finance (see, e.g., [15]), especially in the OU-type stochastic volatility 
model of [6]. An OU-type process is given as the solution of a stochastic 
differential equation of the form 

(1.1) dXt = -aXtdt + dLt 

with L being a Levy process (see, e.g., [45] for a comprehensive introduction) 
and a G ]R. Typically, one is interested mainly in stationary solutions of 
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(1.1). Provided a > and E{ln{\L\ V 1)) < oo, the SDE (1.1) has a unique 
stationary solution given by 



However, in many apphcations the dependence structure exhibited by 
empirical data is found to be not in good accordance with that of OU-type 
processes which have autocorrelation functions of the form e~"'^ for positive 
lags h. In many data sets a more complex and often a (quasi) long memory 
behavior of the autocorrelation function is encountered. OU-type processes 
could be replaced by fractional OU-type processes (see [30] or [31], for in- 
stance) to have long memory effects included in the model. However, in this 
case many desirable properties are lost, and, in particular, fractional OU- 
type processes no longer have jumps. An alternative to obtain long memory 
from OU-type processes and still to have jumps is to add up countably many 
independent OU-type processes, that is. 



with independent identically distributed Levy processes (Lj)jgN and appro- 
priate Oj > 0, Wi > with X^^i^^i = 1- Intuitively we can likewise "add" 
(i.e., integrate) up independent OU-type processes with all parameter val- 
ues a > possible. The resulting processes are called supOU processes and 
have been introduced in [2] where it has also been established that they 
may exhibit long-range dependence. For a comprehensive treatment regard- 
ing the theory and use of univariate supOU processes in finance we refer to 



So far supOU processes have only been considered in the univariate case. 
However, in many applications it is crucial to model several time series with 
a joint model, and so flexible multivariate models are important. Therefore, 
in this paper we introduce and study multivariate supOU processes. Due 
to the appearance of matrices and the related peculiarities our theory is 
not a straightforward extension of the univariate results. Multivariate {d- 
dimensional) OU-type processes (see, e.g., [26] or [46]) are defined as the 
solutions of SDEs of the form 



with L a d-dimensional Levy process and A a d x d matrix. Provided 
i?(ln(||L|| V 1)) < oo and all eigenvalues of A have strictly negative real part, 
we have again a unique stationary solution given by 





oo 



[7]. 



(1.2) 



dXt 



AXt dt + dL; 
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Intuitively our multivariate supOU processes are obtained by "adding up" 
independent OU-type processes with all possible parameters A; that is, we 
consider all d x d matrices A with eigenvalues of strictly negative real parts. 
It turns out later on that the behavior of supOU becomes more tractable 
when we restrict A to come only from a nice subset, like the negative definite 
matrices. 

The remainder of this paper is structured as follows. The next section 
starts with a brief overview of important notation and conventions used in 
the paper and is followed in Section 2.2 by a comprehensive introduction 
into Levy bases and the related integration theory, which will be needed to 
define supOU processes. In Section 3 we first define multivariate supOU pro- 
cesses and provide existence conditions in Section 3.1. Thereafter, we discuss 
the existence of moments and derive the second-order structure. For the fi- 
nite variation case we show important path properties in Section 3.3. Besides 
establishing that we have cadlag paths of bounded variation, we give an ana- 
logue of the stochastic differential equation (1.2) for supOU processes and its 
proof. In particular, this proves a conjecture in [2], which has not yet been 
shown in any nondegenerate set-up. We conclude that section with several 
examples illustrating the behavior and properties of supOU processes and 
showing that they may exhibit long memory. In Section 4 we use our results 
to define positive semi-definite supOU processes and analyze their proper- 
ties. These processes are important for applications like stochastic volatility 
modeling, since they may be used to describe the stochastic dynamics of 
a latent covariance matrix. Finally, this and other possible applications of 
supOU processes are discussed in Section 5. 

2. Background and preliminaries. 

2.1. Notation. We denote the set of real m x n matrices by Mm,ni^)- 
If in = n, we simply write M„(R) and denote the group of invertible n x n 
matrices by GLn(M), the linear subspace of symmetric matrices by Sn, the 
(closed) positive semi-definite cone by §^ and the open positive definite cone 
by (likewise E>~~ are the strictly negative definite matrices, etc.). In 
stands for the nxn identity matrix. The tensor (Kronecker) product of two 
matrices A, B is written as A® B. vec denotes the well-known vectorisation 

2 

operator that maps the nxn matrices to M" by stacking the columns of 
the matrices below one another. For more information regarding the tensor 
product and vec operator we refer to [23], Chapter 4. The spectrum of a 
matrix is denoted by <t(-). Finally, A* is the transpose (adjoint) of a matrix 
A € Mm^„(]R), and Aij stands for the entry of A in the ith row and jth 
column. 

Norms of vectors or matrices are denoted by || • H- If the norm is not 
further specified, then it is understood that we take the Euclidean norm or 
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its induced operator norm, respectively. However, due to the equivalence of 
all norms none of our results really depends on the choice of norms. 

For a complex number z we denote by 3?(z) its real part. Moreover, the 
indicator function of a set A is written 1^. 

A mapping f iV is said to be ^-y^-measurable if it is measurable 

when the a-algebra is used on the domain V, and the c-algebra W is 
used on the range W. The Borel cr-algebras are denoted by and A 

typically stands for the Lebesgue measure which in vector or matrix spaces 
is understood to be defined as the product of the coordinatewise Lebesgue 
measures. 

Throughout we assume that all random variables and processes are de- 
fined on a given complete probability space equipped with an 
appropriate filtration when relevant. 

Furthermore, we employ an intuitive notation with respect to the (stochas- 
tic) integration with matrix-valued integrators referring to any of the stan- 
dard texts (e.g., [43]) for a comprehensive treatment of the theory of stochas- 
tic integration. Let {At)t^^+ in Mm,n(^), (-Bt)jgR+ in Mr^s(M) be cadlag 
and adapted processes and {Lt)t^-g^+ in Mn^r{^) be a semi-martingale. Then 
we denote by Jf^As-dLgBg- the matrix Cj in Mm,s(M) which has ijth 

element Cij^t = Y.k=iYA=i lo '^ik,s~Bij^s- dLki,s- Equivalently such an in- 
tegral can be understood in the sense of [34, 35] by identifying it with 
the integral Jq As- dLg with Aj being for each fixed t the linear operator 
Mn^rO^) Mm^sO^)^^ '"^ AtXBt. Analogous notation is used in the context 
of integrals with respect to random measures. 

Finally, integrals of the form J"^ f{x, y)m{dx, dy) are understood to be 
over the set A'm.x and over B in y. 

2.2. Levy bases. To lay the foundations for the definition of vector-valued 
supOU processes, we give now a summary of Levy bases and the related in- 
tegration theory. In this context recall that a d-dimensional Levy process 
can be understood as an M'^-valued random measure on the real numbers. If 
L = (Lj)jgiR is a d-dimensional Levy process, this measure is simply deter- 
mined by L((a, b]) = L{b) — L{a) for all a, 6 S M, a < 6. 

Define now M~ := {X e Mrf(M) : a{X) C (-oo, 0) + iR} and ^6(M~ x M) 
to be the bounded Borel sets of x M. Note that is obviously a cone, 
but not a convex one (cf. [22], for instance). Moreover, we obviously have 
= {X G Md{R) : a{X) C (-oo, 0] + iR}. 

Definition 2.1. A family A = {A{B) : B G ^^(M^ x M)} of M'^-valued 
random variables is called an R'^-valued Levy basis on x R if: 

(a) the distribution of A{B) is infinitely divisible for all B G =^fe(M^ x M), 
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(b) for any n G N and pairwise disjoint sets Bi, . . . , Bn S ,^b{Mj x R) the 
random variables A(i?i), . . . , A(i?„) are independent and 

(c) for any pairwise disjoint sets Bi G ^i,{M^ x M), i G N, satisfying UnGN £ 
J^b{Mj X M) the series X]5S=i ^{^n) converges almost surely and it holds 
that A(U„eM^n) = E~=iA(i3n). 

In the literature Levy bases are also often called infinitely divisible inde- 
pendently scattered random measures (abbreviated i.d.i.s.r.m.) instead. 

In the following we will only consider Levy bases, which are homogeneous 
(in time) and factorizable (into the effects of one underlying infinitely di- 
visible distribution and a probability distribution on M^); that is, their 
characteristic function is of the form 

(2.1) E{exp{iu*A{B)))=exp{ip{u)U{B)) 

for all ueM'^ and B e ^b(M~(IR) x M). Here 11 = vr x A is the product of a 
probability measure vr on M^(]R) and the Lebesgue measure A on R and 

(2.2) (/?(n) =m*7- ^n*S'u+ / (e*"*^ - 1 - m*xl[o,i](||x||))i/((ix) 

is the cumulant transform of an infinitely divisible distribution on M.'^ with 
Levy-Khintchine triplet (7,S,z^), that is, 7 G M.'^, S G and u is a Levy 
measure — a Borel measure on with i^({0}) = and /jjd(||a;|P A l)i'{dx) < 
00. The quadruple (7,S,z^, vr) determines the distribution of the Levy basis 
completely and is henceforth referred to as the "generating quadruple" (cf. 
[18]). 

The Levy process L defined by 

Lt = A(M^- X {0,t]) and L^t = HM;[ x (-t,0)) for t G M+, 

has characteristic triplet (7, S, v) and is called "the underlying Levy pro- 
cess." 

For more information on M'^-valued Levy bases see [39] and [44]. 
A Levy basis has a Levy-Ito decomposition. 

Theorem 2.2 (Levy-Ito decomposition). Let A be a homogeneous and 
factorisable -valued Levy basis on x M with generating quadruple 
(7,S,i/, vr). Then there exists a modification A of A which is also a Levy ba- 
sis with generating quadruple (7,S,z^, vr) such that there exists an W^-valued 
Levy basis A'^ on x M with generating quadruple (0, S,0,-7r) and an in- 
dependent Poisson random measure fi on (M'^ x x R,^(M'^ x x M)) 
with intensity measure v x ir x X which satisfy 

A(5) =7(7r X X){B) + A^{B) 
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(2.3) + / / x{fi{dx,dA,ds) - dsTT{dA)i^{dx)) 

\\x\\<l Jb 



+ / xfi{dx,dA,ds) 

J\\x\\>i Jb 



for all B G ^b{M^ x M) and all ueQ. 

Provided i||a.||<i ||a;||i^((ix) < oo, it holds that 

A{B)=jo{Tr X X){B) + A^{B)+ [ [ xfi{dx,dA,ds) 

Jr^ Jb 

for all B G ^b{M^ x M) with 



(2.4) 7o:=7- / xu{dx). 

J\\x\\<l 

Furthermore, the integral with respect to fi exists as a Lebesgue integral for 
all uj 

Here an M'^-valued Levy basis A on x M is called a modification of a 
Levy basis A if A{B) = A{B) a.s. for all B £ ^b{Mj x R). For the necessary 
background on the integration with respect to Poisson random measures we 
refer to [24], Section 2.1, and [27], Lemma 12.13. 

Proof of Theorem 2.2. This follows immediately from [39], Theo- 
rem 4.5, because the control measure m is given by m{B) = {\\"f\\ + tr(S) + 
/jgd(l A ||rE|p)i/(dx))(7r x \){B) which is trivially continuous due to the pres- 
ence of the Lebesgue measure. The second part is an immediate consequence, 
as no compensation for the small jumps is needed if J||2,.||<i ||a;||i^(c?x) < oo. 
□ 

From now on we assume without loss of generality that all Levy bases are 
such that they have the Levy-Ito decomposition (2.3). 

In the following we need to define integrals of deterministic functions with 
respect to a Levy basis A. Following [44], for simple functions / : x M ^ 



f{x) = ^ailB,{x) 



i=l 



with m € N, G Mrf(M) and Bi G ^b{M~ x M), and for every B G ^{M~ x 
R) we define the integral 



/ f{x)A{dx) = Y,aiA{Br\Bi). 
-'^ 1=1 
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A =^(A/^ X M)-^(Mrf(M))-measurable function / : A/^ x M ^ M^iR) is said 
to be A-integrable if there exists a sequence of simple functions (/n)nGN such 
that fn^ f Lebesgue almost everywhere, and for all B S SS{M'^ x R) the 
sequence Jg/„(x)A(dx) converges in probability. For A-integrable / we set 

f{x)A{dx) = plim„^oo /„(x)A((ix). As in [44] well definedness of the 
integral is ensured by [48]. 

The following result is a straightforward generalization of [44], Proposi- 
tions 2.6 and 2.7, to M'^-valued Levy bases and follows along the same lines. 

Proposition 2.3. Let A be an W^-valued Levy basis with character- 
istic function of the form (2.1) and f-.M' xR^ Md{R) a Sg{M~ x M)- 
^(Mti{R)) -measurable function. Then f is A-integrable if and only if 



M- 



(2.5) 



f{A,sh 

+ / f{A,s)xil[o,i]{\\fiA,s)x\ 



llo,i]{\\x\\)Mdx) 



ds7r{dA) < oo, 



(2.6) 
(2.7) 



m: 



\\f {A,s)^f {A,s)*\\dsTT{dA) <oo, 



{lA\\f{A,s)xf)v{dx)ds'K{dA) < oo. 



If f is A-integrable, the distribution of Jj^j- Jj^+ f{A, s)A{dA,ds) is 



in- 



finitely divisible with characteristic function 

f{A,s)A{dA,ds) 



E(exp(z.* 
(2.8) 

= exp( [ [ ^{f{A,s)*u)dsTT{dA) 

Vja/^ Jr+ 

and characteristic triplet (7int) ^int, I'int) given by 



7int = I \ fiA,s)-f 



(2.9) 



+ / f{A,s)x{l[o,i]{\\f{As)x\ 



l[o,i]{\\x\\)Mdx)\ds7r{dA) 
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(2.10) Si„t= / / f{A,s)^f{A,srds7r{dA), 

(2.11) Z^int(5)= / / / lB{f{A,s)x)i^{dx)dsTT{dA) V5 G ^(M"^). 

J Af^ Jr 

When the underlying Levy process has finite variation we can do w-wise 
Lebesgue integration; that is, the integral can be obtained as a Lebesgue 
integral for each € 

Proposition 2.4. Let A be an -valued Levy basis with characteristic 
quadruple (7,0,i/,7r) satisfying j^^^^^^_^^\\x\\i'{dx) <oo, and define 70 as in 

(2.4), that is, ip{u) = m*7o + /jg<;(e*"*^ — l)i^((ix). Furthermore, let f-M^ x 
R Md{R) be a ,^{M~ x R)-3g{Md{M.)) -measurable function satisfying 

(2.12) [ [ \\f{A,s)\\dsTr{dA) <oo, 

Jm~ Jr 

(2.13) / / / {1 A \\f{ A, s)x\\)iy{dx)ds 7r{d A) < 00. 
J M- Jr Jr'' 

Then 



(2.14) 



f{A,s)A{dA,ds)= [ [ f{A,shods7r{dA) 
Jmt Jr 



R'i J 



f{A, s)xfi{dx, dA, ds), 



and the right-hand side is a Lebesgue integral for every uj £ Q [conditions 
(2.12) and (2.13) are also necessary for this]. 

Moreover, the distribution of J.,- L /(^, s)A(dA, ds) is infinitely divisible 

d 

with characteristic function 

Elexpliu* I I f{A,s)A{dA,ds) 
\ V Jm^ Jr 

where 



(2.15) 7int,o= / / fiA,s)jods7TidA), 

Jm~ Jr 

(2.16) UintiB) = I [ [ lBifiA,s)x)i^{dx)ds7T{dA) V5 

JM7 Jr JRd- 
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Proof. Follows from the Levy-Ito decomposition and the usual integra- 
tion theory with respect to Poisson random measures (see [27], Lemma 12.13). 
□ 

Remark 2.5. All results of this section remain valid when replacing 
with Mfc(R), € N, or any measurable subset of a finite-dimensional real 
vector space and when considering integration of functions / : (M) x M ^ 
Mm,dO^)- We decided to state all our results with as this set will be 
used mainly in the following and it reduces the notational burden. 

3. Multidimensional supOU processes. In this section we introduce sup- 
OU processes taking values in with d S N and analyze their properties. 
This extends to a multivariate setting the theory of univariate supOU pro- 
cesses as introduced in [2] and studied further, for example, in [19]. 

Intuitively supOU processes are obtained by "adding up" independent 
OU-type processes with different mean reversion coefficient. 

3.1. Definition and existence. We define a d-dimensional supOU process 
as a process of the form (3.4) below. 

Theorem 3.1. Let A be an M.'^ -valued Levy basis on x M with gen- 
erating quadruple (7,S,z/, vr) satisfying 



(3.1) / ln(||a;||)i/((ix) < oo, 

J\\x\\>l 

and assume there exist measurable functions p : — ?• M'^\{0} and k : 
[1, oo) such that 

(3.2) lle^' II < K{A)e-P^'^^' Vs G M+, vr-a/mosi surely, 
and 

(3.3) / ^^7r((iA) < oo. 



JM- P{A) 
Then the process {Xt)t^^ given by 

(3.4) Xt= f f e^^'~''>A{dA,ds) 

is well defined for allt (zM. and stationary. The distribution of Xt is infinitely 
divisible with characteristic triplet {'~fx,^x,'^x) given by 

lx = ( [ fe^^7+ / e^^x(l[o,i](l|e^^x||) 

JM' Jr+ \ JRd 



10 

(3.5) 



O. E. BARNDORFF-NIELSEN AND R. STELZER 



^0,i]i\\x\\)Mdx))ds7ridA), 



(3.6) Sx= / / e'^'J:e^*'dsTT{dA), 
Jmt Jr+ 



(3.7) 



M- JR+ 



lB(e x)i/(dx)(is7r(d^) 



for all Borel sets B C . 



Proof. The stationarity is obvious once the well definedness is shown. 
Using Proposition 2.3 it follows that necessary and sufficient conditions for 
the integral to exist are given by 



(3.8) 



(3.9) 

(3.10) 
First we show (3.10) 



e^«^+ / e^^x(l[o,i](||e^^x|I) 

-l[o,i](||x||))z.((ix) 
/ / ||e^^Se^*^||ds7r(d^) <oo, 

JM^ Jr+ 

(1 A \\e'^^x\\^)u{dx)ds7r{dA) < oo. 



ds Tr(dA) < oo, 



JR+ 



< 



{l^\\e'^'xf)v{dx)dsn{dA) 

{l^K{Afe-'^P^^>\\xf)v{dx)ds^T{dA) 



\n{K{A)\\x\\) + l/2 



+ 



M- J\\x\\>1/k{A) P{^) 

K{Af\\x\\^ 



v{dx)-K{dA) 



M7 J\\x\\<1/k{A) 2p(^) 



v{dx)'K{dA). 



The finiteness of the first integral follows from (3.1), (3.3), k,{A) > 1 and v 
being a Levy measure, which imply 



/ / 

JM^ J\\x\\>1Ik 



>1/k{A) 



ln(K(^)||x||) +1/2 



u{dx)'K{dA) 



< 
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M- J\\x\\>l 



piA) 



u{dx)'K{dA) 



+ 



2|t^l|2 



+ 



+ 



3k(A)2 



-u{dx)7r{dA) 



■K{dA) I v{dx) 

x\\>l 



M- MA) 

1 

A/7 lAA) 



■K{dA) / \\x\Yv{dx) 
Ik-Il<i 



7T{dA) I (ln(||x|l) + l/2)i/(cix) <(X). 

x||>l 



Likewise the finiteness of the second integral is implied by (3.3), k{A) > 1 
and i||2.||<i < oo, as is a Levy measure. 

Next (3.9) follows from (3.3) and 

/ / |le^^Se^*^||ds7r(d^) < ||S|| / / K{Afe-^P^^'>' ds7r{dA) 
Jm~ Jr+ Jm~ Jr+ 

K{Af 



2p{A) 



■n{dA). 



Turning to (3.8) we have from (3.3) that 



e^'-i\\dsTi{dA) < 



M- Jm.+ 



/<A) 
M- P{A) 



K{A)e-P'^^^' dsn{dA) 



■K{dA) < oo. 



Moreover, 



< 



+ 



e^^x(l[_i,i](||e^^rr||)-l[_i,i](||x||))Krfx) 
e"^" x\\v{dx) ds'K{dA) 

-,As 



ds7r{dA) 



M~ JR+ -'||x||<l,||e^''a;||>l 



m: 



x||>l,||e^''x||<l 



e x\\i'{dx) ds7r{dA) 



< 



e'^''x\\^u{dx)ds7r{dA) 



MT JR+ -'||x||<l,||e^''x||>l 
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+ / / / \\x\\K{A)e-P^'^^'iy{dx)ds7r{dA) 

Jm~ Jr+ J\\x\\G{l,eP(A)''/^) 

+ / / / u{dx)dsTr{dA) 

with the finiteness following from (3.1), (3.3) and v being a Levy measure. 

That the distribution of Xt is infinitely divisible and has the stated char- 
acteristic triplet follows now immediately from Proposition 2.3. □ 

Remark 3.2. (i) The necessary and sufficient conditions for the exis- 
tence of the multivariate supOU process X are (3.8)~(3.10). However, as 
they are obviously very intricate to check in concrete situations, it seems to 
be appropriate to replace them by the sufficient conditions (3.1)-(3.3). One 
particular advantage of these conditions is that they involve only integrals 
with respect to either or vr, but not with respect to both. 

(ii) Note also that for d=l the conditions above become the necessary 
and sufficient conditions of [19], as we can then take k,{A) = 1 and p{A) = 
-A. 

(iii) By looking at the Jordan decomposition one can see that pointwise 
there is for any A e a constant k G [1, oo) and a p G (0, — max(5R(cj(A)))] 
such that lle"^*!! < kc"''* for all s G M"*". If A is diagonalizable, it is possible 
to choose p{A) = - max(5R((j(^))) and n{A) = || with [/ G GLrf(C) 
being such that UAU~^ is diagonal. So (3.2) essentially demands that this 
choice has to be done measurably in A (see especially Example 3.5 for a 
concrete example). 

It seems important, especially for applications, to understand how close 
our sufficient existence conditions are to necessity and to give also necessary 
conditions easier checkable than (3.8)-(3.10). To this end we need the con- 
cept of modulus of injectivity (see, e.g., [40], Section B.3). For a Z G M^^W) 
the modulus of injectivity is defined as j{Z) — min||2.||=]^ It is immedi- 

ate to see that 0< j(Z) < \\Z\\ and j{Z) = ||Z~^||~^ for Z G GLd{M.). 

Proposition 3.3. Let A he an W^-valued Levy basis on x R with 
generating quadruple (7,S,z^, vr) and assume there exist measurable functions 
t:M~^ M+\{0} and '&:M~ ^ (0, 1] such that 

(3.11) j{e^') > i?(y4)e-^(^)" Vs G R+,TT-almost surely. 
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Then necessary conditions for the integral (3.4) to exist are 



(3.12) 



/ -^■K{dA) < oo 

J^(A)>e 



h{A)>e '^{^) 

for any e G (0, 1] with y({\\x\\ > l/e}) > and 7r({??(^) > e}) > 0, 
(3.13) f ^^^j-7r{dA) < oo provided > or u{{\\x\\ < 1}) > 



and 

(3.14) / ln(|lx||)i/((ix) <oo. 

J||x||>l 

This shows, in particular, that the logarithmic moment condition on v 
is both necessary and sufficient and that (3.2), (3.3) together with (3.11)- 
(3.13) form a set of sufficient and necessary conditions which are as close as 
one can probably hope for conditions reasonably easy to work with. 

Proof of Proposition 3.3. In the case j(^) > condition (3.13) 
follows immediately from (3.9). 

Equation (3.10) implies that a necessary condition is 

oo> / / / {l^{){Afe-'^^^'^^'\\xf)v{dx)ds^:{dA) 

M^i J\\x\\>l/-9(A) 

+ [ [ ^^f\f u{dxMdA). 

Since '&{A) < 1, the second summand implies (3.13) given i/({||x|| < 1}) > 0. 

Choose now e > such that tt{{'&{A) > e}) > 0. Then the first integral is 
bigger than 

J^{A)>e J\\x\\>l/e '^(^j 

This gives the necessity of (3.12). Moreover, the above integral is again 
greater than 

ln(e) + Indian) + 1/2 



i'{dx)TT{dA). 

^{A)>eJ\\x\\>l/e 

Since /^(^)>£ i||a;||>i/e :f^;)^{dx)TT{dA) is necessarily finite, we have that 

^-^Y^u{dx)7r{dA) > -oo. 

■&(A)>eJ\\x\\>l/e 
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Hence, i||x||>i/e ^-^Tm^''(dx)Tr{dA) < oo. Since /^^^^^^ Tc^^lc?^) > 

by construction, this implies the necessity of (3.14). □ 

Remark 3.4. (i) For d=l we have again recovered the necessary and 
sufficient conditions of [19], as we can then take i!)(A) = 1 and t{A) = — A. 

(ii) Pointwise there is again for any A S a constant € (0, 1] and a 
r € [-min(3fJ(cr(^))),oo) such that i(e^*) > i^e"^^ for ah s S ]R+. For A di- 
agonahzable T(yl) = — min(K((T(^))) and '!9(^) = j{U)j{U~^) can be chosen 
if UAU~^ is diagonal. 

(iii) If has unbounded support, /^(•y^)>£ 7(lj'^('^^) < oo is necessary for 
ah e > 0. If Tr{{'&{A) > e}) = 1 for some e > 0, then /^^- :p^^TT{dA) < oo 
becomes a necessary condition, provided z^({||x|| > 1/e}) > 0. 

In some applications like stochastic volatility modeling, for instance, one 
is particularly interested in the case where the underlying Levy process is 
of finite variation and the supOU process is defined via w-wise integration. 
The following result is proved using Proposition 2.4 together with variations 
of the arguments of the proofs of Theorem 3.1 and Proposition 3.3. 

Proposition 3.5. (i) Let A be an W'- -valued Levy basis on x R 
with generating quadruple (7,0,z/, vr) satisfying 



(3.15) / ln(||x||)i^((ix) < oo and / ||x||z^(da;) < oo 
J||x-||>i •^lla;||<i 

and assume there exist measurable functions p : — > M+\{0} and k : 
[1, oo) such that 

(3.16) ||e^' II < K{A)e~P^'^^' Vs G M+, vr-a/mosi surely, 
and 

(3.17) / '^TT{dA) < oo. 



M 



7 PiA) ' 



Then the process {Xt)t<=R given by 

Xt= I f e^^'~''>A{dA,ds) 

J M7 J-oo 



e^(*-")7ods7r((iA)+ [ [ [ e^^^-''^ xp{dx,dA,ds) 



is well defined as a Lebesgue integral for all t gM. and u £ Q, and X is 
stationary. 
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(ii) // there exist measurable functions r : — ?• M''"\{0} and iD : — > 
(0,1] such that j{e^'^) > 'd{A)e~'^^^^'^ Vs G ,TT-almost surely, then neces- 
sary conditions for the integral (3.18) to exist as a Lebesgue integral are: 

J^A)>e •^(^J 

(3.19) 

for any e G (0, 1] such that i^({||a;|| > 1/e}) > and 7r({i9(yl) > e}) > 0, 

(3.20) [ '^^Tr{dA)<oo provided or u{{\\x\\ < 1}) > 0, 
Jm- r{A) 

(3.21) / ln(||x||)z^((ix) < oo and / ||2;||z^((ix) < oo. 

J||2:||>1 -'bll<l 

Remark 3.6. If (3.3) is satisfied for a Levy basis, then (3.17) is also 
satisfied. 



We shall not develop the general case further, but consider two special 

'TV 
d 



cases which appear to be sufficient for most purposes. We define 



{A G Mrf(M) : A is normal and a{A) C (-oo,0) + 

Proposition 3.7. (i) Assume that tt{M^-) = 1, then (3.2) or (3.16) 
are satisfied with k{A) = 1 and p{A) = — max{^{a{A))) . Moreover, (3.3) or 
(3.17) are implied by 

(3.22) - / .J^. . ^... TT{dA) < oo. 

Aff- max(5?(a(A))) 

Likewise, (3.11) is satisfied with ^{A) = 1 and t(A) = —mm{?R.{a{A))). 
The necessary conditions (3.12), (3.13), (3.19) and (3.20) all become 



( ■ ^^rx/ / .NNN T^{dA) < oo 

Jm:^~ mm{^{a{A))) 



[assuming i/(||x|| > 1) > for (3.12), (3.19)]. 

(ii) Assume that there are a X G N and diagonalizable Ai,. . . , Ak G M^(M) 
such that 7r{{XAi:i = l,...,K;XeR+\{0}}) = 1. Then (3.2) or (3.16) are 
satisfied with k{A) = C for some C G [l,oo) and p{A) = — max(5i(cj(yl))). 
Moreover, (3.3) or (3.17) are implied by 

(3.23) - / — ^— — ^(d^) < oo. 

A/- max(SR(a(A))) 

Likewise, (3.11) is satisfied with t{A) = — mm{^{a{A))) and i?(^) = c for 
some cG (0,1] and the necessary conditions (3.12), (3.13), (3.19) and (3.20) 
all become - Jj^- min(SR(o-(A))) ^(^^) < °° [assuming i^{\\x\\ > 1/c) > for 
(3.12), (3.19)]. 
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In dimension one these are again the weh-known necessary and sufficient 
conditions. Observe also that the eigenvalues are continuous (and hence 
measurable) in A because they are the zeros of the characteristic polynomial. 

Proof of Proposition 3.7. Part (i) follows immediately from the fact 
that all normal matrices are unitarily diagonalizable. 

Likewise, (ii) is a consequence of the above mentioned pointwise bound 
and the fact that this can be turned into a global one because for fixed 
i = 1, . . . , N the matrices {Aylj}AeK+\{o} are all diagonalized by the same 
invertible matrices. □ 

In (i) the mean reversion parameter A of the superimposed OU-type 
processes is restricted to normal matrices and in (ii) to finitely many rays 
{A^i}AeR+\{o}- 

Remark 3.8. (i) Typically one will, in general, not consider normal 
matrices for A as in (i), but only negative definite ones, since this allows one 
to use well-known distributions on the positive definite matrices (see, e.g., 
[21]) for TT. In the case (ii) possible vr can be obtained by using arbitrary 
distributions on along the rays and positive weights summing to one for 
the different rays. 

(ii) Intuitively (3.22) and (3.23) and their necessary counterparts mean 
that vr must not put too much mass on elements with very slow exponential 
decay rates. 

3.2. Finiteness of moments and second-order structure. Before we look 
at the second-order structure, we give conditions ensuring the finiteness of 
moments. 

Theorem 3.9. Let X be a stationary d-dimensional supOU process 
driven by a Levy basis A satisfying the conditions of Theorem 3.1. 



J\\x\\>l 

for r € (0,2], then X has a finite rth moment, that is, E{\\Xt\\^') < oo. 
(ii) // r G (2, oo) and 



(i) If 



(3.24) 




(3.25) 




then X has a finite rth moment, that is, EdlXtW^) < oo. 
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(iii) Necessary conditions for X to have a finite rth moment are 

(3.26) / ||x|rz^((ij;) <oo 

J\\x\\>l 

in general and 

(3.27) / ^AdA) < oo 

for any e such that z^({||x|| > l/e}) > and 7t{{'&{A) > e}) > 0. 

In connection with the above results observe that the underlying Levy 
process L has an rth moment, that is, E'dlLi H*") < oo, for r € if and only 
if 



/ ||x||^i''L((fx) is finite. 

J||x||>l 



Proof of Theorem 3.9. Using [45], Corollary 25.8, we have to show 
i|bll>i lkir^x(^^) < oo to establish (i) and (ii). Now, 



\\x\\''' vxidx) 
x||>i 



M- JO 



e x\ 



'Mi,oo){\\e^"x\\)u{dx) dsTT{dA) 



< / / K{AYe-''P^^^'\\x\\''l(^-^^^){K{A)e-P^^^'\\x\\)u{dx)dsTr{dA) 
Jm^ Jo Jwi 

ln{K(A)\\x\\)/p{A) 

K{AYe-''P^^>\\x\\'' ds u{dx)TT{dA) 



K{AY\\xf[^ 1 

J\\x\\>1/k{A) 

K(^)'^||x|r-1 

M- J\\x\\>i/k{A) rp{A) 



^^^T^ f 1 - ^TT^^l u{dx)TT{dA) 
rp{A) \ K{Ay\\x\\'' J ^ ^ ^ ^ 

u{dx)TT{dA). 



That /^j- i||x||>i/K(A) rp{A) ^{dx)iT{dA) < oo has already been shown in the 
proof of Theorem 3.1. 
Moreover, we obtain 



k{AY\\x\ 

M~ J\\x\\>1/k{A) 



-v{dx)'K{dA) 
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Jm- J\\x\\>1 P(^) 

/4NrV2|| ||rV2 

"^^^ .idx)nidA). 

M^; J\\x\\<i P{^) 

Hence, (i) and (ii) follow, since i/ is a Levy measure, using also (3.3) for (i). 
Regarding the proof of (iii), analogous arguments give that 

' i^{dx)TT{dA) < oo 

M- J\\x\\>i/-&{A) rT[A) 

is a necessary condition. An inspection of the proof of Proposition 3.3 shows 
that 



{l/T{A))i^{dx)TT{dA) < oo 

J\\x\\>l/-&{A) 

is already necessary for X to exist. Hence, 

oo > / / , u{dx)-K{dA) 

Jm~ J\\x\\>i/^A) rT{A) 

> [ [ "^^^^[[f ,.{dx)7r{dA) 

J'9(A)>eJ\\x\\>l/e rT{A) 

is necessary for a finite rth moment of X. This implies (iii), since 
/||x||>i ll'^ir^(^-^) finite if and only if ll-^ir^('^-^) ^ arbitrary 

OO. □ 



Remark 3.10. In the set-up of Proposition 3.7(i) [and analogously in 
(ii)] 

(3.28) -/ —^———TT{dA)<oo, I ||xf i/((ix) <oo 

Jm- max(K(cj(A))) y||^||>i 

imply (3.22) and (3.25) [resp. (3.24)]. 

Likewise, - /^,^- min(jR(<7(A))) ^(^^) < oo, provided u{\\x\\ > 1) > or iy{\\x\\ > 
1/c) > 0, respectively, and ||a;||''i/((ix) < oo become necessary condi- 

tions for X to exist and to have a finite rth moment. 



In applications these results have important implications for modeling. If 
one wants to have certain moments finite and certain moments infinite, for 
example, because this is what observed data strongly suggests, one has to 
use a driving Levy basis A having exactly the same moments finite. 
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Moreover, knowledge of the moments allows for the estimation of the 
model based on empirical observations by using the general method of mo- 
ments (GMM) estimation procedure, for example. Often the inference on 
parameters is based on the second-order moment structure as, for instance, 
in the estimation of multivariate OU-type stochastic volatility models in [41] 
or the one of univariate supOU type models in [7]. To provide the foundations 
for such work and due to the importance in applications of understanding the 
temporal dependence structure, we now calculate the second-order moment 
structure. 



Theorem 3.11. Let X be a stationary d- dimensional supOU process 
driven by a Levy basis A satisfying the conditions of Theorem 3.1 and assume 
additionally that f^a \\x\\'^ij{dx) < oo. Then E{\\Xq\\'^) < oo and we have 



(3.29) E{Xo) = - [ A~^(j+ [ xi^{dx))Tr{dA), 

JM7 V J\x\>l J 



xx*v{dx) ]-K{dA) 



(3.30) var(Xo) = - / {s^{A)r^ i^+i I 

coY{Xh,Xo) = - [ e^'^(^(^))-i(s+ / xx*u{dx))ir{dA) 

(3.31) 



for h € . 



with £^{A) : Md(M) Md(M),X ^ AX + XA* . 
Moreover, it holds that 

(3.32) lim cov(Xft,Xo) =0. 



Proof. The finiteness of the second moments follows from Theorem 
3.9. Using the formulae of Theorem 3.1 and [45], Example 25.12, we obtain 

E{Xo)='yx+ / xux{dx)= / / e^'*(7-F / xiy{dx)] ds7r{dA). 

J\\x\\>l Jm~ Jr+ \ -'||a:||>l / 

Noting that ■^A~^e'^^ = e"^^, integrating over s gives (3.29). 
Likewise we get 

var(Xo) = Sx + / xx*ux{dx) 

= [ / e^"(s + / xx*u{dx)]e'^*'ds7r{dA) 
Jm- Jr+ \ JR'i ) 

which implies (3.30) by integrating over s. 
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Finally, 



cov{Xh,Xo) = cov( [ _ r e^^^-''^A{dA,ds), [ _ T e-^'AidA,ds)) 

V^A/^ J~oo JrI^ J~od J 

[ f e^^''-'^A{dA,ds), [ f e-^'A{dA,ds) 
f e-^^fs+ / xx*u{dx)\e-'^*'d^T:{dA) 



(3.33) 



e^'^(=c/(^))-M S+ / xx*u{dx)]TT{dA) 



since A is a Levy basis, and hence the random measures A on x (0, h] 
and on x (— oo,0] are independent. 
From (3.33) one obtains 



Ah 



m: 



f e~^^fs+ / xx*u{dx)\e-^''' d^^{dA) 



< 



< 



J-oo 



xx*v{dx) 



m: 



2p{A) 



7r{dA) 



S+ / xx*v{dx) 



< oo. 



Therefore lim/i_!.oo = for all A G and dominated convergence es- 
tablish (3.32). □ 



3.3. "SDE representation" and some important path properties. In this 
section we show for a supOU process X a representation which generalises 
the SDE that governs OU-type processes, and we derive important path 
properties of X. The "SDE representation" — identity (3.38) below — has 
been conjectured in the univariate case in [2], where neither a proof nor 
conditions for its validity have been given. Below we are able to show these 
results for finite variation Levy bases, which are naturally appearing in appli- 
cations like stochastic volatility modeling. The properties which we establish 
are especially important in the context of integration, since they imply that, 
if X is the integrator, then pathwise Lebesgue integration can be carried out, 
and, when X is the integrand, the theory of stochastic integrals of cadlag 
processes with respect to semimartingales (see [43], for instance) respectively 
the L^-theory of, for example, [38] applies. Likewise, the integrated process 
is of importance in certain applications (see Section 5.2). 
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Below the filtration {J^t)teM generated by A is defined by J^t being the 
(j-algebra generated by the set of random variables {h.{B) : B € ^{MJ x 
(-00, t])} for i gM. 

Theorem 3.12. Let X he a supOU process as in Proposition 3.5. Then: 

(i) Xt{uj) is ^(M) X measurable as a function 0/ t G M and w G 
and adapted to the filtration {,^t)teR generated by A. 

(ii) // 



< 00, 



(3.34) / K{A)Tr{dA) 

the paths of X are locall-^ uniformly bounded in t for every a; G fi. 
Furthermore, X^ = Jq Xg ds exists for all t G M"*" and 

X+= [ f A~^e^^'~''^A{dA,ds)- [ f A-^e-'^'KidA.ds) 
A^^A{dA,ds). 



M-Jo 



(iii) Provided that 

f (ll^ll Vl)K(yl) 

(3.36) / ^" '' , ^ \ (dA) < 00 

Jm- p{^) 



and 

(3.37) / \\A\\K{A)Tr{dA) < 00 
it holds that 

(3.38) Xt = Xo+ [ Z^du + Lt, 

Jo 

where L is the underlying Levy process and 

(3.39) Zu= [ r Ae^^''-'^A{dA,ds) 

for all u (zM with the integral existing u-wise. 

Moreover, the paths of X are cddldg and of finite variation on compacts. 

Proof, (i) is immediate from the definition of Xt as a Lebesgue integral 
and the measurability properties of the integrand e^^^~^^xl^+{t — s) which 
as a function of t,s,A,x is =^(M x M x x ]R'^)-^(R'^)-measurable. 



22 O. E. BARNDORFF-NIELSEN AND R. STELZER 

(ii) We first show local uniform boundedness of X. Choose arbitrary 
Ti, r2 G M with Ti<T2. Then 

fTuT2{A,s,x):= sup ||e^(*-")xlK+(i-s)|| 

iG[Ti,T2] 

for all A G M^(M), s G M and x eR'^ and 



sup \\Xt\\< / / fTi,T2{A,S,Jo)dsTT{dA) 



m: 



/ti,T2 (^) ■5, x)fi{dx, dA, ds). 



Therefore we only have to show the w-wise existence and finiteness of the 
integral on the right-hand side. This is, however, an immediate consequence 
of the above upper bound, (3.34), Proposition 2.4 and arguments as in the 
proof of Theorem 3.1 noting that 

{1 A K{A)\\x\\)u{dx) ds7r{dA) 



<{T2-Ti)( I [ K{A)\\xMdx)Tr{dA) 
\Jm^ ■jM\<i 

lu{dx)TT{dA) 

-'||2;||>l 

Turning to the existence follows immediately from the local bound- 
edness. Noting that we have actually proved the local boundedness of 

e^^*-'^jo\\dsTr{dA)+ [ f f We^^^-'') x\\^i{dx,dA,ds) 

above, we can use Fubini to obtain 

X+= I f f e'^'^'-''^-fodsduTr{dA) 

^11 I I e'^^'~''''xdsfi{dx,dA,du) 

= I f A-\^^^-^ho\Uoyu)dunm 

+ / / _ r A-'e^^'-^^x\l^^^^^^dsfi{dx,dA,du), 
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which estabhshes (3.35) by straightforward calculations. 

(iii) Using similar calculations as before, the existence of as an w-wise 
integral follows from Proposition 2.4 and (3.36). Similarly to (ii) one sees 
that under (3.37) Z is locally uniformly bounded in u. Hence, one can use 
Fubini to obtain 

ct r r rt ft 



I 

Jo 



Zudu = [ [ I I Ae^'^''-'^xdu^i{dx,dA,ds) 

Jm.'' JaI^ J~oo Jo\/s 

+ I f f Ae^^'^-'hQdudsT:{dA) 

Jm~ J -oo JoVs 



rt 



x\u=(ovs)Kdx,dA,ds) 



+ I [ e^^""'HolL(ovs) ds^dA) =Xt-Xo-Lt 
Jm- J-oo 

which establishes (3.38). That X has cadlag paths of finite variation is now 
an immediate consequence of this integral representation. □ 

Remark 3.13. (i) Condition (3.34) is always true if tt is concentrated 
on the normal matrices or on finitely many rays and hence especially in 
dimension d = 1. Moreover, it could be replaced by the weaker but rather 
impracticable condition that /[^^^ A 1 is integrable with respect 

to TT X A X z/ and that, for any fixed x, f[Ti,T2]i^j ^^x) is integrable with 
respect to vr x A (cf. [29], Proposition 2.1, for a very related result whose 
proof is similar in spirit to ours, but uses a series representation instead of 
the Levy-Ito decomposition). 

(ii) Intuitively (3.37) means that vr does not place too much mass on the 
elements of with high norm and thus very fast exponential decay rates. 

If vr is concentrated on the normal matrices or finitely many diagonalizable 
rays, then (3.36) and (3.37) become 

(3.40) - / vr(dA) < oo and f \\A\\7r{dA) < oo. 

JAq maxK(o-(A)) J m- 

In particular, the second condition simply means that vr has a finite first 
moment. 

If vr is concentrated on then we have ||^|| = — min(cr(A)), and (3.36) 
becomes 

/ . X f (min(cr(A)) A -1) , , 

3.41 / ^ \ T^{dA) < oo, 

max(a(^)) 

so it can be seen as a condition on the spread between the different ex- 
ponential decay rates measured by the eigenvalues. It is easy to see that 
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in dimension d=l, it is equivalent to J^^{ — l/A)-K{dA) < oo, which is part 
of the necessary and sufficient conditions for the existence of the supOU 
process. 

(iii) It is very easy to construct examples when our sufficient conditions for 
X to exist as an w-wise integral are satisfied, but neither (3.36) nor (3.37) 

for Z above. Take, for example, vr concentrated on Vn = (^~f^ with 

T^ivn) = 6/(7r^n^). In such a case we unfortunately do not know whether 
Z exists because our previously employed techniques seem, at best, to give 
a necessary condition of the type (3.36) involving j{A) A 1, and hence the 
necessary conditions for Z to exist would be implied by the sufficient ones 
for X. Therefore, we also refrained from giving necessary conditions in this 
section. 

3.4. Examples and long-range dependence. Like in the univariate case, 
the expression (3.31) does not imply that we necessarily have an exponential 
decay of the autocovariance function and thus a short memory process. On 
the contrary we can easily obtain a long memory process, as the following 
examples exhibit. Note that this illustrates that (3.32) is not obvious and 
indeed requires a detailed proof as above. 

Apart from showing that multivariate supOU processes may exhibit long- 
range dependence, the purpose of this section is to analyze some concrete 
examples and their properties. 

Regarding long-range dependence, there is, unfortunately, basically no 
general theory developed in the multivariate case so far. Below we mean by 
long-range dependence (or long memory) simply that at least one element 
of the autocovariance function decays asymptotically like h~°' for the lag 
h going to infinity and for some a G (0;1). Intuitively this should clearly 
be a case when one may appropriately speak of long-range dependence. 
Establishing a general theory for multivariate long-range dependence seems 
to be very important, but is beyond the scope of this paper. 

Example 3.1. Let A be a d-dimensional Levy basis with generating 
quadruple (7, S, z/, vr) with satisfying ||x|pi/((ix) < 00 and vr being given 
as the distribution of RB with a diagonalizable B G and R being a real 
r(a, /3)-distributed random variable with a > l,/3 € M"*"\{0}. Hence, R has 
probability density /(r) = j^^'^"~^e~^'"l]g+ (r), and from 
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max(K(f7(S)))r(Q) amax(5R(a(5))) 

we conclude that (3.28) holds. Hence, the process Xt = Jjyj- e^^^~^^ A{dA, 

ds) exists, is stationary and has finite second moments. Similar calculations 
imply that a > 1 is also necessary for Xt to exist. 

For the autocovariance function at positive lags h we find 



cov{Xh,Xo 



m: 



xx*v[dx\ 



UBU- 



with ^:Md(]R) Md{M),X ^ BX + XB*. Let now U G GLd(C) and Ai, A2, 
. . . , Ad e (—CO, 0) + zM be such that 

/Ai ••• 0\ 
A2 '•• i 

: ■•. ■•. 

VO ••• XdJ 

Then, from t^~^e~^^ dt = r{z)k~^ for aU z,k G (0,oo) + iM, where the 
power is defined via the principal branch of the complex logarithm (see [1], 
page 255), we obtain that 



u 



exp 



/Ai 



A2 



o\ W 

h 



T{a-l)U 



\0 
/Ai 
A2 



V 



■•. 

XdJ 
o\ \ 

h 



r'^~'^drU~^ 



■•. 

XdJ 



u- 



J 



Vo • 

= r(a-l)(/3/rf-S/i)i-". 

Above the (1 — a)th power of a matrix is understood to be defined via 
spectral calculus as usual. 
Hence, 

cov(X;„Xo) = ^(/3/d-5/i)i-°^-ifs+ / xx*u{dx)], 
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and thus we have a polynomiany decaying autocovariance function. For a € 
(1,2) we obviously get long memory. 

Another question is whether we have the nice path properties of The- 
orem 3.12. Hence, assume additionally that i||2,||<i ||x||z^((ix) < oo. In our 
example condition (3.34) is trivially satisfied and so the paths of X are lo- 
cally uniformly bounded in t. Regarding condition (3.40) the second part is 
equivalent to 

/•oo 

||S|| / r'^ e- ^"^ dr <oo, 
Jo 



which is always true, as any Gamma distribution has a finite mean. Denoting 
the density of the r(a, /3)-distribution by /a,/3(?"), one obtains for the first 
part 



\\B 



max{R{a{B))) 



which is obviously finite. Hence, the conditions of Theorem 3.12(iii) are 
satisfied and thus the paths are cadlag and of finite variation, and (3.38) is 
valid. 



Example 3.2. The previous example has an immediate extension to the 
case when tt is concentrated on several rays instead of a single one as above. 
Assume we have wi, . . . ^Wm G [0,1] with Y17!=i'^i — 1 ™d diagonalizable 
i?i , . . . , Bm G , and define vTj to be the probability measure of the random 
variable RiBi with i?j being V[ai,(3i) distributed with Oj > 1, /3i G M'''\{0}. If 
V is as above and tt = Yll=i '^i'^i-, we get for the multivariate supOU process 
X 

coy{Xh, ^o) = - ^ ~ B^hf~"'^i^^ + j^^ xx*u{dx)^ 

with : Mrf(M) ^ Md(M) , X ^ BiX + XB* . 

Assuming now J||^.||<i ||2;||z^((ix) < oo, it is likewise straightforward to see 
that conditions (3.34) and (3.40) are satisfied. Hence, the paths of X are 
locally uniformly bounded in t, cadlag and of finite variation, and (3.38) is 
valid. 



Example 3.3. A similar result can be obtained if we restrict the mean 
reversion parameter A to the strictly negative definite matrices and 
define vr as a probability distribution on the proper convex cone ST~ as 



MULTIVARIATE SUPOU PROCESSES 



27 



follows. Let denote the intersection of the unit sphere in Sd with S^~, 
let a:S^~ (l,oo), /3:S^~ (0,oo) be measurable mappings and w a 
probability distribution on such that 

(3.42) - / ^^^1 w{dv) < oo. 



Now define vr via 



/3(t;)"(^) 



Vs^- r(Q(t;)) 

for any Borel set B G MJ (M) . Then vr is a probability distribution concen- 
trated on 

Using this vr in the above set-up means that the mean reversion parameter 
is no longer necessarily restricted to finitely many rays. Moreover, similar 
calculations to the ones in Example 3.1 give 



^ -7r{dA)= f . . ^^""^ . .. w{dv) <oo. 



m: 



max(JR{a{A))) — a{v) max{a{v)) 



Hence, (3.28) holds and the process Xt = fj,r- fl^^^^'^ '^^A{dA,ds) exists, 

d 

is stationary and has finite second moments. Likewise we get for the auto- 
covariance function 



S -|- / xx*v{dx 



cov{Xh,Xo) 



with y{v):Md{R)^MdiM.),X^vX + Xv*. 

Turning to the path properties, assume now that i|jx.||<i ll^ll'^(f^2;) < oo. 
Again condition (3.34) is trivially satisfied and so the paths of X are locally 
uniformly bounded in t. Regarding condition (3.40) the second part becomes 

[ \\A\\7r{dA)=[ I rfa{v),i3[v){r)drw{dv) = I ^^^w{dv), 
Jm- Js-- Jo Js-- P[v) 

and for the first part one obtains 

A/- maxK(cr(A)) ^ ^ 

(3.43) = - / / ]——fa(v),i3{v){r)drw{dv) 

1 



' / / ■.■, fa{v),i3{v){r)drw{dv). 

I max(cj(f)j ^ ' ^ ' 
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The first summand is finite due to (3.42) and the second one is finite if the 
integral — L — {1/ 'niayi{a{v)))'w{dv) is finite. Hence, provided 

d 

-w{av) < oo and / w[av) < oo, 



max{a{v)) J^-- 

the conditions of Theorem 3.12(iii) are satisfied, and thus the paths are 
cadlag and of finite variation and (3.38) is valid. 

Based on this we can easily give an example where we know that the 
supOU process exists due to Proposition 3.5, but the conditions of Theorem 
3.12(iii) are not satisfied. Assume w is a discrete distribution concentrated 
on the points 



n G N, 

and that 'w{vn) = -r^, a{vn) = 2 and (3{vn) = n~^. Then we have that 




but 



/ 

J A 



12 °° 



A\\Tr{dA) = n ^ = oo, 



n=l 

and hence condition (3.37) is not satisfied. Observe that this means that the 
probability measure vr we have constructed does not have a first moment, 
although it is defined via a polar representation where the radial parts are 
all univariate Gamma distributions. Moreover, it is easy to see that (3.43) 
is finite and so Z exists and thus it is only the local uniform boundedness 
our sufficient conditions fail to provide when trying to show (3.38). Showing 
that (3.38) is indeed not valid seems to be a very delicate issue, as already 
remarked. 

Example 3.4. Let A be now a two-dimensional Levy basis with gener- 
ating quadruple (7, S,z^, vr) with v satisfying jj^a ||a;|pi/((ix) < 00. We restrict 
the mean reversion parameter A to B^~, the 2x2 diagonal matrices with 
strictly negative entries on the diagonal. Hence, vr is a measure on '^2~ 1 
which can be identified with (R )^, and we assume that vr has Lebesgue 
density 

■K{dai,da2) 

- ^"'^2' , -ai)"i-i(-a2)"^-^e^^'^i+/^^'^n(K-p(ai,a2)daida2 



r(ai)r(Q2) 
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with ai,Q2 > 1 and /3i,/?2 > 0. So the diagonal elements are independent, 
and their absolute values follow Gamma distributions. We obtain 

^ -Tr{dA) 



max(3?(cj(A))) 

Dai 002 



fco foo 1 oai 0012 

Jq Jq m.m{ai,a2)T{ai)T{a2) 



oai poo 002 

Jo r(ai) Jq T{a2) 

POO pOLX i-OO Q012 

Jo r(Qi) Jq r(Q2) 

Hence, (3.28) holds, and the process Xt = fj^,^~ e'^^*~**-*A((iA, ds) exists, 
is stationary and has finite second moments. 

Let us now consider the individual components Xi^t, -^2,t of Xf. Denote by 
Pi : ^ R, (xi, X2)* I— T- xi the projection onto the first coordinate and define 
an M- valued Levy basis Ai on M x M via Ai{dai,ds) = P{A{P^^{dai),ds) 
and a Levy measure ui on M via ui{dxi) = u(P^^{dxi)). Then Ai has char- 
acteristic quadruple (71, Sn, z^i, tti) with vri having Lebesgue density 

r(ai)) ^ ^ 

and 



Xi,t= / / e'^i(*-^)Ai((iai,ds). 
Jr — J-00 

For the autocovariance function of the first component we get 

cov(Xi,fe, Xi,o) = ' + f Sn + / xjuiidxi)), heR+. 

2(ai - 1) \ Jr ) 

An analogous result holds for the second component X2^t and we have long 
memory in both components provided 01,02 G (1)2). 

The importance of this example is, however, that we can model the sta- 
tionary distributions of X\ and X2, that is, the margins of the stationary 
distribution of Xf, very explicitly by specifying the margins of i/, that is, 
v\ and V2- From [2], Theorem 3.1, Corollary 3.1, and [19], Remark 2.2, we 
know that it is exactly all nondegenerate self-decomposable distributions 
on M which arise as the stationary distributions of the components. More- 
over, these authors provide formulae to calculate v\ (or V2) if one wants 
to obtain a given stationary distribution for the component (alternatively 
[3], Lemma 5.1, or the refinement [42], Theorem 4.9, can be used). Hence, 
one can specify a two-dimensional supOU process with prescribed stationary 
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distributions of the components by calculating the required vi and 1/2 and 
choosing v accordingly. The easiest way to get a possible v is by specify- 
ing v{dxi,dx2) = vi{dxi) X 5q{x2) + 5q{xi) x 7/2(^x2) with 5q denoting the 
Dirac distribution with unit mass at zero. In this case the components of 
X are independent. An easy way to get an appropriate v and allowing for 
dependence is to combine ui and 1^2 using a Levy copula (see [4, 28]). 

Likewise it is again interesting to look at the path properties of Theorem 
3.12. Assuming again i||j.||<i |lx||z^(d2;) < 00, condition (3.34) is trivially sat- 
isfied, and so the paths of X are locally uniformly bounded in t. Regarding 
the second part of condition (3.40) we have that 

/ \\A\\T:{dA)<oo 

is equivalent to 

/ max(ai, a2)/ai,/3i {ai)fa2,p^{a2) dai da2 

J(R+)2 

< / aifa^,i3^{ai)dai+ / 02/0,3,^2(02) da2 < 00, 
which is always true. Turning to (3.41), it is implied by 



/ 



max(ai,a2) , . / ^J ^ 

—f r/ai,/3i(,aij/a2,/32l02j «Oi (202 

m.m.(ai,a2) 

< / "^"^ /ai,/3i(ai)/a2,fe(°2)c^Ql ^Q2 



+ / -^/ai,/3i(ai)/a2,/32(a2)d«lC?a2 < OO, 

J(R+)2 0,2 

which is easily seen to be always true. Hence, the conditions of Theorem 
3.12(iii) are satisfied, and thus the paths are cadlag and of finite variation 
and (3.38) is valid. 

Obviously this example has a straightforward extension to general dimen- 
sion d. 



Example 3.5. So far we have only studied cases where we could use 
Proposition 3.7 and did especially never have to bother with k{A) in the 
conditions of Theorem 3.1. 

In this example we will present a case where the behavior of k.{A) is crucial 
and where we show how k and p can be specified in a measurable way. We 
define the following sets: 

= {X € M^(M) -.X is diagonal; all diagonal elements are strictly negative. 
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pairwise distinct and ordered such that xa < xjj VI < i < j < d}; 
= {X € GL(i(]R) : the first nonzero element in each column is 1}; 

= {SDS-^ ■.Seyd,D£ &-}. 

If A = SDS~^ is in , the matrix D consists of the eigenvalues of A, 
and the columns of S are the eigenvectors of A. In principle there are many 
possible S and D if we only demand A = SDS~^. However, if we restrict 
ourselves to S* € G , then S,D are unique, as elementary linear 

algebra shows. This means that the map 

m-.^dX &~ {S, D) ^ SDS-^ 

is bijective (and obviously continuous). We denote by lUl"^ = (©,3) the in- 
verse mapping. Since computing eigenvectors and eigenvalues are measurable 
procedures as are the orderings and normalizations involved in obtaining the 
diagonal matrix in and the eigenvector matrix in S^d, all these mappings 
are measurable. Note also that ^d-, Borel sets. 

Defining K-.Ji' [l,oo),A ^ ||e(^)|| ||(6(^))-i || , p{A) = 

— max(5R((T(^))) gives, therefore, measurable mappings on satisfying 
lle^'^ll < K{A)e~P^'^^'^ . Using these definitions for k and p one could now spec- 
ify probability distributions vr on and check whether condition (3.3) is 
satisfied, and the associated supOU process therefore exists. 

However, in concrete situations it seems easier to specify a Borel proba- 
bility measure T^y^^r^- on .y^d x ^d define vr as its image under 9Jt, i.e. 

7^(3) = vTy^^^- {Tl~^{B)) for all Borel sets B. Assume T^y^^r^- = T^y^ ^ t^^- 
is the product of two probability measures tt on yd and ir^- on S'^ . Then 
we have 

/ ^^71T^('^^) < ^ / \\Sf\\S~'^fTTy^{dS)<oo and 
Jai- P{A) Jy^ 

J^- max(3fe(cj(iJ))) -^d 

That \\S\\'^\\S~^\\'^TTy^{dS) can be finite or infinite depending on the 
choice of ny^ is exhibited by the following example. Let ny^ be a discrete 
measure concentrated on the points 

-Sn = ^ ^ 1 ) ^^'^ ^" ^'^2 i^n) = Can~°' Vn G N 
with a > 1 and Ca = ^ / Yl'^=i^~"' ■ Then 

s„- = ( > » 

-n 1 
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Using the equivalence of all norms we get that 

r. oo 

/ |lS'|p||S'~"'"|p7r^2(dS') < cxo <^=^ Ca/ n^Pn<oo <;=^ a > 5. 

Returning to the general example with vr given via vr^^ x vr^- and turning 
to path properties, we assume again i||^||<i ||a;||z^((ix) < oo. In this finite 
variation case the existence conditions (3.17) become 

/ ||5|| IIS""^ ||7rv'j(c?5') < oo and — / , Tr^^- (dD) < oo. 

Jy, ' J 92 max(K(a(Z)))) ^ > 

Furthermore, condition (3.34) is always satisfied when the existence condi- 
tions are satisfied and so the paths of X are locally uniformly bounded in 
t. Straightforward arguments show that the conditions of Theorem 3.12(iii) 
are satisfied, and thus the paths are cadlag and of finite variation and (3.38) 
is valid if 

|2|| c-l||2„ ' ^ 11-^1 



5f |lS"^f7r,y',(dS) <oo, -/ __t^,l,r^^^-, T^oj-{dD)<oo 



^'-""-^ J^-max(K(a(Z))))"^. 
and 



/ ||D||7r^-(dZ)) < oo. 

J917 



By polarly decomposing vr^- into a measure on the unit sphere in the di- 
agonal matrices and a radial part, the long memory specifications of the 
foregoing examples have straightforward extensions to this set-up. 

4. Positive semi-definite supOU processes. Based on the previous sec- 
tion we now consider supOU processes which are positive semi-definite at 
all times. The importance of such processes is that they can be used to de- 
scribe the random evolution of a latent covariance matrix over time and, 
hence, they can be used in multivariate models for heteroskedastic data, for 
example, the stochastic volatility model of [9]. 

Let us briefiy recall that a d x d positive semi-definite OU-type process 
(see [8]) is defined as the unique cadlag solution of the SDE 

dj:t = iAj:t + j:tA*)dt + dLt, Sogs+, 

with A S Mrf(M) and L being a dx d matrix subordinator (see [5]), that is, a 
Levy process in with Lt — Lg^ Vs, t G M"*", s < t. If £^(ln(max(||Li ||, 1))) < 
00 and max.(JR[a{A))) < 0, the above SDE has the unique stationary solution 
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That the hnear operators — ^ of the form Z ^ AZ + ZA* with some A € 
Mrf(]R) are the ones to be used for positive semi-definite OU-type processes 
has been estabhshed in [42]. 

As just recalled, one has to restrict the driving Levy process to matrix 
subordinators in order to obtain OU-type processes taking values in the posi- 
tive semi-definite matrices. Below we need to impose a comparable condition 
on the Levy basis to get positive semi-definite supOU processes. Note that 
for a d X d matrix- valued Levy basis A we denote by vec(A) the M*^^ -valued 
Levy basis given by vec{A){B) =vec{A{B)) for all Borel sets B. Moreover, 
observe that tr{XY*) [with X,Y & M^(M) and tr denoting the usual trace 
functional] defines a scalar product on Mrf(M) and that the vec operator is 
a Hilbert space isometry between Mrf(M) equipped with this scalar product 
and M"^^ with the usual Euclidean scalar product. 

Positive semi-definite supOU processes are defined as processes of the 
form (4.4) below, which is the analogue of (3.4). 

Theorem 4.1. Let A be an Sd-valued Levy basis on x M with gen- 
erating quadruple (7,0, z^, vr) with := ^ — J^^^^f^^^xv{dx) G and v being 

a Levy measure on satisfying u{^d\^~d) — 0; 



(4.1) / ln(||x||)zv((ix) < 00 and I \\x\\u{dx)<oo. 
J||x||>i ■^lkll<i 

Moreover, assume there exist measurable functions P-M^ — )• M+\{0} and 
n'-M^ [l>oo) such that 

(4.2) lle^'^ II < K{A)e-p'^^^' Vs G R+,TT-almost surely, 
and 

(4.3) / ^^7r(dA) < 00. 



P{A) 



Then the process (Sj)j(=ir given by 

Et= / r e^^'-'^A{dA,ds)e^*'^'-'^ 

(4.4) = / f e^^'-'hoe^'^'~'^dsT:{dA) 



J -00 

+ [ [ f e^'^'-'^xe'^*'^'-'^fi{dx,dA,ds) 
JSd Jm~ J -00 

is well defined as a Lebesgue integral for all t gM. and u G ft and S is 
stationary. 



34 



O. E. BARNDORFF-NIELSEN AND R. STELZER 



Moreover, 

(4.5) vec(St)= / f e(^®^'*+^'*®^)(*-^)vec(A)(dA,ds), 

St G for all t and the distribution of Sj is infinitely divisible with 
characteristic function 

S(exp(itr(7xSt)))=exp(^ar(n7E,o) + ^ (e* l)z.E(dx)^ , ueSa, 
where 

(4.6) 7s,o= / / e^^7oe^*^(isvr(d^), 

/■ roo r 

(4.7) i^e(S)= / / / lB(e^"xe^*")i/(dx)ds7r((i^) 

/or all Borel sets B C S^. 

Proof. The equivalence of (4.5) and (4.4) follows from standard results 
on the vectorization operator and the tensor product (see [23]). 

Next we note that e(^®^d+/d®A)(t-s) = gA^^A ^^laX. ||e^(8)e^|| = ||e^||2 
(using the operator norm associated with the Euclidean norm). Hence, all as- 
sertions except St € S^^ for all t E ]R+ follow immediately from Propositions 
2.4 and 3.5. 

However, G §^ for all t G is now immediate, since the integral exists 
cj-wise, e^^Xe"^** G S+ Vvl G Mrf(M), X G § s G M and S+ is a closed convex 
cone. □ 

Remark 4.2. (i) As in Proposition 3.7, k{A) can be replaced by 1 and 
p{A) by -max(3?(cT(^))) in (4.3) [and also in (4.14) and (4.15) below] pro- 
vided TT is concentrated on the normal matrices or finitely many diagonaliz- 
able rays. 

(ii) Throughout this section we refrain from stating necessary conditions, 
as they can be immediately inferred from the foregoing sections and the 
arguments presented for the sufficient conditions. 

Most importantly in the context of stochastic volatility models, which 
involve stochastic integrals with S as integrand. Theorem 3.12 also has an 
analogue for positive semi-definite supOU processes. 



Theorem 4.3. Let S be the positive semi- definite supOU process of 
Theorem 4-1- Then: 
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(i) is =^(]R) X ^ measurable as a function o/t G M and u; € and 
adapted to the filtration {,!^t)teR generated by A. 

(ii) // 

(4.8) / K{AfTr{dA) < oo, 

the paths of T, are locall'^ uniformly bounded in t for every a; G O . 
Furthermore, Tif = Jq ds exists for all t G M"*" and 

S+= / r (A(A))-i(e^(*^")A(dAt^s)e"^*(*-^)) 

(4.9) - I _ f (A(A))-i(e-^"A(d^,ds)e-^*") 

I {K{A)Y^K{dA,ds) 

M- Jo 

with A{A) ■.Sd-^Sd,X ^ AX + XA* . 

(iii) Provided that 

f (ll^ll V1)k(A)2 
4.10 - / ^" " ' n{dA) < oo 

Jm- p{^) 



and 

(4.11) / \\A\\K{A)^Tr{dA) <oo 
it holds that 

(4.12) St = So+ / Zudu + Lt, 

Jo 

where L is the underlying matrix subordinator and 

(Ae^(^-")A((iA,ds)e^*("-") 

MT J-oo 

(4.13) 

+ e^("-'^)A(dA, ds)e^*^''-'^A*) 

for all u (zM with the integral existing uo-wise. 

Moreover, the paths ofTi are cddldg and of finite variation on compacts. 

Formula (4.9) is of particular interest in connection with stochastic volatil- 
ity modeling, as in this case the integrated volatility is a quantity of 
fundamental importance (see Section 5.2). 

Finally, we consider the existence of moments and the second-order struc- 
ture which follow immediately from Theorems 3.9 and 3.11. 
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Proposition 4.4. Let T, be a stationary S'^ -valued supOU process driven 
by a Levy basis A satisfying the conditions of Theorem 4-1- 
(i) If 



(4.14) / ||x||''z^((ix) < oo 

J\\x\\>l 

for r € (0, 1], then T, has a finite rth moment, that is, i?(||S(||'') < oo. 

(ii) //r€(l,oo) and 

(4.15) [ ||2;|ri/(dx) < oo, / ^^"1^], Tr(dA) < oo, 

J\\x\\>l Jm- P{A) 

then S has a finite rth moment, that is, E^WTitW^) < oo. 

(iii) // the conditions given in (ii) are satisfied for r = 2, then the second- 
order structure of S is given by 



i?(So) = - / A{Ar'(jo+ I 



xv{dx) I -K^dA) 



var(vec(So)) = - / (i^(A))"M/ vec{x)vec{xyi^{dx)]Tr{dA), 



cov(vec(E;,),vec(So)) = - / e^^^^''+i'i^^)h(.^{A))-^ 



Yec{x)YGc{x)*v{dx) \ TT{dA) 
d / 

with A{A) : Md(M) ^Md{M.),X^ AX + XA* and £/{A) : M^2 (M) Mrf2 
X^{A(^Id + Id0 A)X + X{A* /rf + /rf O A*). 

Examples 3.1-3.5 can all be immediately adapted to the positive semi- 
definite set-up. More examples in connection with stochastic volatility mod- 
eling can be found in [9]. 

5. Areas of applications. In this section we discuss possible applications 
for our model and the relevance of our results for them. Some of these 
applications are already developed further in other work. 

5.1. Time series modeling. In many areas of applications (e.g., telecom- 
munication, hydrology, economics, finance) one is confronted with time series 
exhibiting a long memory behavior (see [16], for instance) or at least a decay 
of the autocovariance appearing to be a polynomial decay rather than the ex- 
ponential as typically encountered in models of Markovian nature. Moreover, 
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often the relevant data series are multidimensional, and multivariate models 
are needed in order to understand and adequately model the dependence 
effects of the observed data. For irregularly-spaced or high-frequency data 
as well as when intending to look at the data at more than one frequency, it 
is often advisable not to use discrete-time models (like AR(FI)MA, see, e.g., 
[12]), but continuous-time models. Such models can exhibit both continuous 
and discontinuous sample paths. Among the ones with continuous sample 
paths are continuous-time counterparts of ARFIMA, like FICARMA (see 
[13, 31, 47]), which also exhibit long memory. Often it is, however, appro- 
priate to use models with discontinuous sample paths. In such situations it 
appears adequate to use multivariate supOU processes. It should be noted 
that the individual autocovariances of multivariate supOU processes do not 
necessarily have to decay monotonically like in the univariate case, but may 
exhibit damped sinusoidal-like and comparable behavior due to the involved 
matrix exponentials. Hence, one can reproduce second-order moment be- 
havior which in the univariate case calls for the use of CARMA(p, g')-based 
models. 

As for all such models it is a delicate issue to estimate multivariate supOU 
processes from observed data, and this poses many challenging questions. 
Our calculation of the second-order moment structure clearly opens the door 
for (general) method-of-moments-based techniques. Let us illustrate this in 
a simple example. 

Consider the set-up of Example 3.1 with B restricted to have only real 
eigenvalues. Note that in this set-up the parameters cannot be identified, 
because replacing /? by /3c and B by cB leaves the law of the Levy basis 
invariant for any c > 0. So we assume (3 = 1 without loss of generality. Recall 
that 

acov(/i) :=cov(Xfe,Xo) = - ^"^'^~^^^^ [ xx*u{dx) 

a - 1 V JR'i 

with .^:Mrf(M) Md(M),X BX + XB*. Assuming that S + J^^ xx*v{dx) 
is invertible we can define 

Th = acov(/i) acov(O)-^ = {h - Bhy-"" 

which has eigenvalues of the form f{h) = (1 — A/i)^~" with A € (— oo,0). So 
we can obtain an estimator d by calculating the empirical autocovariance 
function from data and fitting / to the maximum (or any other) eigenvalue 
of aSov(/i)aSov(0)~^ for h = nA with n € N and A > being the distance 
between subsequent observations (assumed to be constant) by nonlinear 
least squares, for instance. Thereafter, we can get an estimator for B hj B = 
Id — (ciSov(/i)acov(0)~^)"~^. Now it is straightforward to also get estimators 
for J + Ju^n^-^xi'{dx) and S + J^axx*i'{dx) from the empirical mean and 
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variance. One thus obtains the first two moments of the underlying Levy 
process. If we restrict the ahowed Levy bases such that the first two moments 
of the underlying Levy process identify the parameters 7, S, v, we have thus 
obtained a procedure to estimate all parameters of our model. 

To prove consistency and asymptotic normality (or other asymptotic dis- 
tributions in the case of true long memory) of the estimators one obviously 
needs to understand the (highly non-Markovian) dependence structure of 
our processes and establish mixing conditions or appropriate substitutes. 
We hope to address this issue in future work. 

In a set-up like in Example 3.4 the estimation becomes much easier be- 
cause all parameters except the ones describing the dependence between the 
components can be inferred from the univariate marginal distributions and 
the univariate moment structure. In particular, the parameters /3j and can 
be estimated from the autocovariances of the individual one-dimensional se- 
ries; hence, one does not have to compute eigenvalues or matrix powers as 
above. 

Finally, it should be noted that in Theorem 3.1 or formula (2.8) we have 
given formulae for the characteristic function of multivariate supOU pro- 
cesses. Hence, one can also use estimation techniques based on the empirical 
characteristic function. This is obviously best done in cases where the driv- 
ing Levy basis is chosen such that the integrals in Theorem 3.1 or formula 
(2.8) can be calculated analytically. Moreover, these formulae can be used 
to calculate higher order cumulants for the methods of moments based esti- 
mation. 

5.2. Finance and econometrics. Levy-based stochastic volatility models 
are very successfully applied in both financial mathematics and financial 
econometrics, since they capture many of the stylized facts (nonconstant, 
stochastic volatility exhibiting jumps, heavy tails, volatility clustering, lever- 
age effect,...; see, e.g., [15, 20]) of financial returns very well. One model 
often employed is the Ornstein-Uhlenbeck-type stochastic volatility model 
introduced in [6]. The use of supOU processes as the volatility process, that 
is, the process modelling the instantaneous (co)variance, allows one to in- 
troduce also long memory, which is another important stylized fact, but not 
covered by most models, into the model. Let us illustrate this in a simple 
set-up. Let S be a positive semi-definite supOU process as introduced in 
Section 4 and satisfying the conditions of Theorem 4.3(ii) and (iii). Then 
the log-returns of d financial assets (stocks or currencies, for instance) are 
given by 



(5.1) 
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with /U,/? G W^, initial log prices Yq independent of A, piS^ — ?• M'^ a linear 
operator and L being the underlying Levy process. It should be noted that 
in this model jumps in the price and the volatility always occur together 
which is reasonable for financial data (see [25]). 

An extension of the above model has been investigated in-depth in [9] 
where it is, in particular, shown that long memory in S causes long-range 
dependence in Y and explicit formulae for the moment structure of the 
(squared) returns are obtained. This allows an in-depth econometric analy- 
sis and estimation of the supOU stochastic volatility model comparable to 
what has been done in [41] for the multivariate OU-type stochastic volatility 
model, where it was shown that that model can be estimated and fits well 
to observed data, both from the stock and foreign exchange markets. 

In financial mathematics one is often interested in calculating prices of 
derivatives from a given model and in determining the parameters from 
option prices observed on the markets, referred to as calibration. In most 
reasonably realistic models, unlike the Black-Scholes model, one cannot ob- 
tain closed-form formulae for the prices of derivatives. However, calculating 
prices via Monte Carlo simulations is typically too time consuming. When- 
ever possible, better techniques to calculate derivative prices, that is, con- 
ditional expectations of future payoffs, are called for. One technique which 
proved to be very adequate in many situations is the calculation of the prices 
by inverting the Laplace transform (see [14, 17, 37]). For the multivariate 
OU model this technique is successfully applied in [36]. 

Also, in the above given model (5.1) one can calculate the conditional 
Fourier transform of the future prices. Since some complex values will arise, 
one has to be careful with the scalar products. We use the scalar product 
{xi,X2) = X2X1 on M'^ and {Xi,X2) = tr(X2Xi) on Sd denoting by x* the 
Hermitian and by tr the trace of a matrix. For a linear operator * also 
denotes the adjoint operator in the following. Moreover, we assume that 
£^(exp(itr(A(i3)*n)) = exp{ip/^{u)Il{B)) for all u G with ip being the cu- 
mulant transform of the underlying Levy process and 11 = vr x A. Actually, 
E{ex.p{itr{A{B)*u)) exists for all u G Md(M) -|- zS^ and if can be extended 
to this domain as well. Then it follows by similar arguments as in [36, 41] 
that 




(5.2) 
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p u 



ds7r{dA) 



for all ueR'^ and t € (0, oo). Here A{A)~* := {A{A)~^)* is the linear opera- 
tor on Mrf(R) given hy X A*X + XA. If one restates the above formula by 
representing Lp\ in terms of the Levy-Khintchine triplet, one can calculate 
some of the integrals with respect to ds in the drift and Brownian covariance 
matrix part explicitely, since f^^e'^''^^~'^\A{A)~* X)e^^^~^^ ds = X , for ex- 
ample. However, since this results in rather lengthy formulae, especially in 
the part coming from the Levy measure, we refrain from giving further de- 
tails. Note that we only condition on Yq, since Sq is highly non-Markovian 
and, hence, not informative regarding future values of S, and that equation 
(4.9) is essential to obtain (5.2). 

However, it seems to be an important question what to condition upon in 
such a non-Markovian setting and to the best of our knowledge this issue has 
not been addressed so far. In (5.2) we basically assume that we only know 
the current price. One could also assume that one knows all historic prices 
and, hence, the historic values of S, since they are given by the continuous 
quadratic variation of the prices. Unfortunately, it seems extremely hard 
to understand what happens if one conditions on all these historic prices. 
Another point of view would be to say that W and A resemble the informa- 
tion arriving at the markets and that market participants observe all this 
information precisely. In that case it is appropriate to condition upon the 
(T-algebra % generated by A up to time zero and Yq (which now is only 
assumed to be independent of future values of A and W) and one obtains 



£;(e^^**"|%) 



: exp< I 



(5.3) 



{Yo+p,t)*u 

+ trf/ r A{Ay\e^^'-''^AidA,ds)e 

\Jm- J-oo 

- e-^'A{dA, ds)e-^'') (ufi* + ^-uu* 



+ 



A{A)-* ( n/3* + -uu' 



p u 



dsir{dA) 



for ah ueR'^ and t G (0, oo). 
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It is clear that under appropriate technical conditions the conditional 
Laplace transform given Yq or ^q, respectively, exists in a neighborhood of 
zero and (5.2) or (5.3) can be extended to hold on this neighborhood. Then 
one can use Laplace transform techniques to calculate prices of financial 
derivatives. Like for the standard OU-type model in [36] specifications are 
called for under which some of the integrals can be calculated explicitly, 
since otherwise the numerical integration takes too long to make pricing 
and especially calibration feasible in reasonable time. Observe that in (5.3) 
one would set 

Zt:= f f A(A)'i(e^(*-^)A((iA,ds)e^*(*-^) -e-^^A(dA,ds)e-^*") 

and determine Zt also by calibration to option prices. Obviously, one can do 
this only for derivatives with a fixed maturity t £ M"*", unless one increases 
the number of parameters one calibrates. 

5.3. Multivariate supCAR(M^). Ornstein-Uhlenbeck-type processes are 
a special case of the so-called (multivariate) continuous time autoregres- 
sive moving-average (CARMA) processes (see [10, 11, 33]). As the contin- 
uous time analogue of ARM A processes, CARMA processes, are a funda- 
mental class of processes for time series modeling in continuous time. A 
(i-dimensional supCAR(p) process Y can be defined by 

Yt = {IdA---,0) [ f e^^'-'\0,...,0,IdfA{dA,ds), 

dp 

where A is an M'^-valued Levy basis on x M with generating quadruple 
(7, S,z^, vr) with vr concentrated on the matrices in of the form 

/ Id ••• \ 



(5.4) 



: ■•. ■•. ■•. 

Id 

\-Ap -Ap_i ••• -A2 -AiJ 



with appropriate dxd matrices Ai . Clearly, "/l^^ e^^*"'^) (0, . . . , 0, Id)^A{dA, 
ds)" is for A fixed a CAR(p) process (note that this is in contrast to [32] 
who define supCARMA processes differently), so it is appropriate to call 
this process supCAR(p). Obviously many properties for Y follow from our 
results immediately, since it is basically given by the first d-coordinates of 
a high-dimensional supOU process, and this definition gives a possibility to 
extend CAR processes allowing for long memory and jumps. 
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Using our techniques one can define supCARMA(p, q) process with q <p 
using an R*^- valued Levy basis on x M^p^d x IK with generating quadruple 
(7, E,i/, tt). To obtain proper supCARMA processes one demands 7r{£/ x 
Mdp^d) = 1) denoting the set of matrices of the form (5.4) by s^, and one 
sets 



Adapting and extending our arguments one can easily obtain results for this 
class of processes. For the interpretation as CARMA processes note that the 
moving average coefficients have to be calculated from the d x d blocks of 
B by inverting the formulae given in [33], Theorem 3.12. 

6. Conclusion. In this paper we introduced multivariate supOU pro- 
cesses and obtained various important properties of them. Furthermore, 
some areas of application have been outlined and we are currently consid- 
ering their use in stochastic volatility modelling beginning in [9]. However, 
there are still many important issues concerning the supOU processes them- 
selves which we hope to address in future work. Of particular interest is, for 
example, the development of good estimators for supOU models and to show 
properties like consistency and asymptotic normality for them. This is re- 
lated to understanding better the dependence structure of supOU processes, 
which are clearly not Markovian. 

Likewise, we have shown that supOU processes allow to model long mem- 
ory effects (in a specific sense). However, a detailed theory of multivariate 
long-range dependence needs to be developed. 
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